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Abstract 

We present results of a lattice QCD application of a coordinate space 
renormalization scheme for the extraction of renormalization constants 
for flavour non-singlet bilinear quark operators. The method consists in 
the analysis of the small-distance behaviour of correlation functions in 
Euclidean space and has several theoretical and practical advantages, 
in particular: it is gauge invariant, easy to implement and has relatively 
low computational cost. The values of renormalization constants in 
the X-space scheme can be converted to the MS scheme via 4-loop 
continuum perturbative formulae. Our results for Nf = 2 maximally 
twisted mass fermions with tree-level Symanzik improved gauge action 
are compared to the ones from the RI-MOM scheme and show full 
agreement with this method. 

PACS numbers: 11.15.Ha, 12.38.Gc 



1 Introduction 

For many physical quantities to be computed in Lattice QCD, renormaliza- 
tion is an essential ingredient. Therefore, it is extremely important to have, 
ideally, several different methods available for the determination of renormal- 
ization constants (RCs) that are moreover non-perturbative. In this paper, 
our main aim is to show that the extraction of RCs from the behaviour of 
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correlation functions in coordinate space (termed X-space method below) is 
feasible and, besides providing the RCs in a gauge invariant manner, has 
certain advantages. Therefore, the X-space method has the potential as an 
useful alternative to other widely used methods, although, as we will show 
in this paper, further improvements are still necessary to allow for a precise 
computation of RCs with fully controlled systematic uncertainties. 

One way to extract RCs for Lattice QCD is Lattice Perturbation The- 
ory (LPT), see Ref. [1] for a comprehensive review. While with LPT it is 
difficult to go beyond 1-loop ordeiQ it is still very useful, if not necessary to 
complement non-perturbative computations of RCs. One of developments 
in LPT is a tadpole resummation method proposed in Refs. [7J, [E], which 
has been further explored in Refs. [9], [10] to encompass effects of various 
gluonic and fermionic discretizations. Tadpole resummation amounts to a 
redefinition of the coupling constant. A method to go beyond the 1-loop (or 
sometimes 2-loop) approximation is the framework of Numerical Stochastic 
Perturbation Theory (NSPT) [11], [12], |13| . This method, originates from 
the Stochastic Quantization theory of Parisi and Wu [2] , [15] , which can be 
used to formulate Stochastic Perturbation Theory, which is then applied to 
stochastically evaluate perturbatively accessible quantities, such as RCs or 
the running coupling constant by numerical methods. This makes it possible 
to go to high loop orders in LPT. The method has been successfully applied 
in the context of renormalization, see e.g. Refs. [16], [IT], [15] . 

Of course, eventually a non-perturbative approach to determine RCs is 
desirable. The most natural approach is the use of Ward identities. Such 
an approach allows to obtain a very good precision for a selected number 
of RCs, see e.g. Refs. [E], [20], [H], [22], for an introduction see Ref. [23]. 
However, the application of Ward identities is rather limited and cannot be 
used for a number of important renormalization constants. Therefore, in 
order to compute a broad range of RCs one has to resort to different but 
still non-perturbative methods. 

The two most widely used non-perturbative methods are the Regularization- 
Independent Momentum Subtraction (RI-MOM) and the Schrodinger Func- 
tional (SF) renormalization schemes. For an extensive review of their prac- 
tical usage, see Refs. [S^, [25]. The RI-MOM method [26], also called the 
Rome-Southampton method, is similar in the definitions of renormalization 
conditions to continuum perturbation theory. It consists in the lattice evalu- 
ation of quark propagators of the relevant operators between external quark 
and gluon states and requires hence to fix the gauge. Since the original 
method has been proposed, it was subject to many improvements of its 
reliability, including the perturbative computation of the relevant 0(a 2 g 2 ) 

An example of a 2-loop computation of RCs for clover fermions is given in Refs. [2J, 
[3]. 2-loop LPT computation of quark mass renormalization is presented in Ref. [4]. 2- 
loop relation between the bare lattice coupling and the MS coupling in pure SU(N) gauge 
theories can be found in Refs. [5], [B]. 
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discretization effects [27] , [28] and Goldstone pole subtraction [29] , [30J . The 
RI-MOM approach has been applied in several lattice setups. Since in the 
end we want to compare our results for the RCs with the ones computed 
in the RI-MOM scheme using configurations generated by the European 
Twisted Mass Collaboration (ETMC) |31j . [32], we mention here the work 
related to this setup, i.e. Refs. [33], [22], [33]. Other lattice setups included 
e.g. Nf = 2 + 1 + 1 MTM fermions on Iwasaki gauge action [5S] and Nf = 2 
clover fermions on Wilson plaquette action [36J. For completeness, we note in 
passing the approach of the spectral projectors method, see Ref. [37], which 
can be used to compute the scale independent ratio Zp/Z$ and which has 
already been applied to twisted mass fermions. 

The other widely used non-perturbative method, the Schrodinger func- 
tional [38] , [39] is defined in a four-dimensional Euclidean space with Dirich- 
let boundary conditions in time direction. Combined with finite-size tech- 
niques [40] , it can be used to formulate a SF renormalization scheme |41j . 
with the spatial lattice size used to set the (inverse of) renormalization scale. 
In the SF scheme, the scale dependence and the continuum limit are sepa- 
rated by determining the continuum limit of the lattice step scaling function 
first, before the scale dependence is considered. The SF scheme is formu- 
lated in a completely gauge invariant manner, however, carrying through 
the corresponding renormalization programme bears some complexity. The 
method has been successfully applied e.g. in Refs. [12], [13], for an introduc- 
tion see e.g. Refs. [33], [15]. Concerning maximally twisted mass fermions 
as we will address here, the SF scheme formalism has been developed in 
Refs. [16], [37], leading to the introduction of chirally rotated boundary 
conditions. For first applications, however in the quenched theory only, we 
refer to Refs. 08], 09], [50]. 

Our present paper concerns yet another approach to non-perturbative 
renormalization. The X-space method was suggested in Ref. [JT] and the 
results of its first implementation were presented in Refs. [52], [53], j53J. It 
consists in imposing a renormalization condition directly on gauge invariant 
correlation functions at small Euclidean distances, leading thus to the ad- 
vantage of a gauge invariant determination of RCs. A further advantage is 
the absence of contact terms. In principle, this makes the X-space method 
also suitable for renormalization of four-fermion operators, relevant e.g. for 
the computations of weak matrix elements, although in this paper we limit 
ourselves to RCs of bilinear quark operators. For the X-space method, the 
natural limitation is the necessity of existence of a renormalization window, 
namely the regime of Euclidean distances large enough that uncontrollable 
discretization effects can be avoided, but at the same time small enough to 
make contact with (continuum) perturbation theory (i.e. much smaller than 
the inverse of Aqcd)- The purpose of this paper is to explore the X-space 
method in practice and in particular to see, whether in current simulations, 
taking Nf = 2 LQCD setup as an example, the above mentioned window 
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exists and whether the method is feasible. By comparing to results for RCs 
using the RI-MOM scheme for the Nf = 2 maximally twisted mass setup, 
we can also check the validity of the approach. We will also make an at- 
tempt to determine the systematic errors appearing in the computation of 
RCs with the X-space method, which will allow us in the end to address its 
competitiveness with the methods described above. 

The paper is organized as follows. In Sec. [5J we present the theoretical 
foundations of the method. Details of its practical implementation are given 
in Sec. whereas our final results are presented in Sec. [H Sec. [S] concludes 
and discusses the prospects for broader application of the method. 

2 Theoretical foundations 

In this section, we briefly describe theoretical aspects of the X-space method 
and our lattice setup. We start by summarizing twisted mass Lattice QCD, 
then we introduce the correlation functions whose numerical analysis will 
be presented in Sees. HI Next, we provide the definition of the X-space 
renormalization scheme and discuss some of its properties. Finally, we quote 
some results from continuum perturbation theory which are needed for the 
conversion of our results to the MS scheme and to the evolution of the 
renormalization constant to the scale of 2 GeV. 

2.1 Twisted mass Lattice QCD 

We consider correlation functions calculated from gauge link configurations 
generated with Nf = 2 dynamical quarks by ETMC |31| . [32]. The action 
S[tp,tjj,U], S[ip,'tp,U] = Sg[U] + Sf{i/j,i/j,U] used to generate these con- 
figurations is given by the tree-level Symanzik improved gauge action [55J 

S G [U] = ^(b £ Tr(l-P lxl (x; //,!/))+&! £ Tr(l - P lx2 (*;^, „))) , 

(1) 

with (5 = 6/g 2 , g - bare coupling, &o = § an d b\ = P lxl , P lx2 denote 
the plaquette and the rectangular Wilson loop, respectively; the Wilson 
twisted mass fermion action is given in the twisted basis by [56], [57] . [58], 

m 

SF[ip,ip,U] = a A ^2i/j(x)(D w + i/x 9 7 5 r 3 )V'(x), (2) 

X 

with 

D w = \ (7m(V m + V* ) - aV^V^) + mo- (3) 

As usual, in the above formulae a denotes the lattice spacing, whereas V M 
and are the discretized gauge covariant forward and backward deriva- 
tives. The bare untwisted and twisted fermion masses are denoted by mo 
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Un-twisted case Twisted case 





= Z A A fl 




= Z v A fl 


(Vr)„ 


= z v v ll 


(Vr)» 


= z A v fl 


(Pr) 


= ZpP 


(Pr) 


= ZpP 


(Sr) 


= z s s 


(Sr) 


= z s s 



Table 1: Summary of our convention of naming the correlation functions 
and their corresponding renormalization constants. 

and n q , respectively. ip{x) is a two-component vector in flavour space build 
of two spinors u(x) and d(x). The matrix r 3 in Eq. ([2]) acts in this flavour 
space. Such formulation allows for an automatic 0{a) improvement of phys- 
ical observables, provided the hopping parameter k = (8 + 2amo) _1 is tuned 
to maximal twist by setting it to its critical value, at which the PCAC quark 
mass vanishes [56], [60], [20J. 

2.2 Correlation functions 

Correlation functions considered in this work are constructed from flavour 
non-singlet bilinear quark operators and are of the form 

C rr (X) = (Or(X)O r (0)), (4) 

where 

O r (X)=i,(X)Ti;(X) (5) 

with Or = {S, P, Vp, A^} for T = {1, 75, 7 M , T/xTs}- X denotes the Euclidean 
vector: X = (xx,x<2, £3, £4) = (x, y, z, t). 

In order to avoid confusion, it is important to state explicitly our conven- 
tions of denoting the correlation functions and their corresponding RCs. We 
label RCs according to the un-twisted Wilson case notation. Therefore, by 
Zy and Z A we denote the RCs of the local vector and axial- vector currents 
of the standard Wilson action, although in our case they renormalize in the 
twisted basis the local axial-vector and polar-vector currents, respectively 
(see Tab. HI). 

2.3 X-space renormalization scheme 

The renormalization constants are defined non-perturbatively by the fol- 
lowing condition [5l] imposed in position space and in the chiral limit: H 

hm(0*(X)0*(O))| M2 = (Or(X )Or(0)}^ t (6) 

2 The renormalization point will always be denoted by Xo = (xo, yo, zo, to). By Xq we 
mean XE = (xo) 2 + (y ) 2 + (zq) 2 + (to) 2 . When the meaning of Xo is unambiguous, we 
will also use it to denote Xo = \J (xo) 2 + (yo) 2 + (zo) 2 + (to) 2 - 



5 



where the renormalized operator is 



O£(X,X ) = Z£(Xl{x ,yo,2o,k),a)Or{X), (7) 

and Xq is the renormalization point. The vector (xq, yo, in, to) is simply- 
given by -%-{xq, yo, Zo, to) and denotes the direction of the renormalization 

pointd. Xq must satisfy the condition a <C lo « ^qcd so * na ^ * ne dis- 
cretization effects and the connection to continuum perturbation theory are 
under control. 

The condition of Eq. ([6]) should be understood in the following way: For 
every finite value of the lattice spacing a, we impose the condition 



free,massless 
cont ) 



(zX(xl(x ,y ,z ,i ),a)) (O r (X )O T (0))(a) = (O r (X )O r (0)) 

from which the RC at this value of a, Z^(Xq, (xq, yo, Bo, to), a), can be cal- 
culated. The continuum limit of Eq. ([6]) is then trivially satisfied. Such 
procedure then guarantees that we recover the correctly normalized quanti- 
ties of the continuum theory. In order to illustrate this last statement, let us 
consider the isovector current V?. In the continuum theory it is a conserved 
quantity defined by the symmetry transformations 

Si/, = % -a a T a ^, 

6$ = -i> % -a a T a . (9) 

On the lattice, two discretizations are commonly used, which have the same 
naive continuum limit, namely, (see e.g. Ref. [19J) 



1 
2 

1 / 



V£{X) = -j(X)T a ^(X), (10) 

v; 



r «(X) = - (${X) ( 7 „ - 1) U,(I)r>(I + M )+ 

+ $(X + M ) (7m + 1) Ul(X)r a ^{X)) , (11) 

but differ in that is conserved by the lattice version of symmetry trans- 
formations Eq. (|9|), whereas V£ is not. We can relate these two currents by 

a finite normalization factor denoted by Zy s defined as 

v; = z™v«, (12) 



3 In a practical application of this method, we will use an empirical criterion (so called 
"democratic" points) to select such direction (xa,yo, zo,io) that is least affected by lattice 
artefacts and average over few neighbouring directions in order to obtain a stronger signal. 
Details of our procedure will be described in Sec. [3] A systematic analysis of this aspect 
will be presented elsewhere. 
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where we choose the MS scheme only for definiteness. In the X-scheme, one 
has a similar condition 

ZZVZ = Z*V». (13) 

In this case, however, even the conserved current receives a renormalization 
since the renormalization condition Eq. ([6]) is incompatible with the sym- 
metry given by Eq. (J9j) and thus the non-renormalization theorem does not 
hold. Note that this fact already appears in continuum perturbation theory 
and is hence no lattice specific effect. However, since the two discretizations 
given by Eqs. (|10p and (|11|) agree in the continuum limit, the conversion 
to the MS scheme through a perturbative conversion factor is the same for 

2.4 Perturbation theory results 
2.4.1 Correlation functions 

At weak coupling, the functional form of the continuum counterparts of the 
correlation functions Crr(X) introduced by Eq. (jl]), which are denoted by 
lip in the continuum theory, is known up to fourth order in the coupling 
constant [ST] , namely 



oo 

U P (X) = U S (X) = ( 1 + Y, C%%) , (14) 



n=l 



7T 4 (X 2 ) 



(15) 



with 

oo oo 

C v = l + £c^a» D v = Y,Dl~a n s , (16) 

n=l n=0 

where the values of the perturbative coefficients C~[, and CY for n < 4 



are given in Ref. [61]. The equalities of the pseudoscalar and scalar, as well 
as the vector and axial correlation functions, follow from the assumption of 
quarks being massless and of working in the flavour-charged sector in which 
the disconnected graphs are absent. In the above expressions: 

ajfi = l/X) 
a s = , (17) 

7T 



where d s is defined in the MS scheme [61] which is related to the MS scheme 
by a shift of the renormalization scale: 

P- = Mms = 2e " yE f J 'MS ~ L12 / x ms = L12 A- (18) 
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The relation between MS quantities and their MS counterparts is 

IT^pf; p, = l/X) = U m (X; = le^jj), (19) 

and 

o a (A) = a s {j2 = 2e^ B p). (20) 
a s is defined as a s /vr = g 2 /(2n) 2 with 

/4 a ^ = o./3(«.) = -£fta*- a (21) 

" i>0 

and the beta function in the MS scheme is defined by 

oo 

P(a s ) = -Y,M + \ (22) 

n=0 

with the 4-loop /3 n coefficients computed in Ref. [62J. 

2.4.2 Perturbative conversion from the X-space scheme to the 
MS-scheme 

The natural scale for the transition between the MS scheme and the X-space 
scheme is 

= -4- (23) 

This is equivalent to the transition between the MS scheme and the X-space 
scheme at a scale 

Ao = 4e" 2 ^. (24) 
For the correlation functions, we can write 

which yields the desired ratios needed for the conversion between the differ- 
ent schemes: 

Z^(Ao) _ Z^Oio) _ 



n MS (Ao) 



(26) 



I free 



z x (x ) z x (x ) \J n(x ) 

In Ref. |61j . these ratios have been computed perturbatively: 

^4 = i+E^(A r (27) 



n=l 

and 

^ s (Ao) 



Zf(Xo) 

7 x (Y \ = 1 + 2^S n a s ( f i ) 



(28) 



with known numerical values of the coefficients 5^ and 8^ for n < 4. 
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K 


L A 


x T 


(if i 


Label 


a [fm] 


L [fm] 


3.9 


0.160856 


24 3 


x 48 


0.0040 


B40.24 


0.0790(26) 


1.9 










0.0064 


B64.24 














0.0085 


B85.24 














0.0150 


B150.24 






4.05 


0.157010 


32 3 


x 64 


0.0030 


C30.32 


0.0630(20) 


2.0 










0.0060 


C60.32 














0.0080 


C80.32 






4.20 


0.154073 


48 3 


x 96 


0.0020 


D20.48 


0.05142(83) 


2.4 


4.35 


0.151740 


32 3 


x 64 


0.00175 


E17.32 


0.0420(17) 


1.3 



Table 2: Ensembles and their parameters used for the calculations. Re- 
sults for the lattice spacing are taken from Ref. [US] (ensembles B,C,D) and 
Ref. [66] (ensemble E). 

2.4.3 Perturbative evolution of the renormalization constants 

Given a renormalization constant at a given scale \x we can obtain pertur- 
batively its value at another scale // through the renormalization group 
equation (staying solely in the MS-scheme) 

c r Ifi) 

where the perturbative expansion of the function c r /IS (//) for V = {S, P}, re- 
lated to the quark mass anomalous dimension, was computed in Refs. |63| . 
|64j . For r = {V, ^4}, the Ward identities for the axial and vector currents 
imply that their renormalization constants in the MS scheme are scale in- 
dependent, which is equivalent to saying that their anomalous dimensions 
vanish identically. 

3 Analysis strategy and example 
3.1 Lattice setup 

We apply the X-space method of extracting RCs for Nf = 2 dynamical 
ensembles, generated by ETMC [31], [32], [65], using the tree-level Symanzik 
improved gauge action and 2 flavours of twisted mass fermions. 

We consider four values of the inverse gauge coupling /?, corresponding 
to lattice spacing values of between ~ 0.04 and 0.08 fm. The parameters of 
our ensembles are summarized in Tab. [2j 
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We perform our investigation for four types of correlation functions: 



Css(X) = (S(X)S(O)}, Cpp(X) = (P(X)P(O)), 
C VV (X) = ^(V^(X)^(O)), C AA (X) = £<A M (X)A M (0)>., 

which allows us to extract the four renormalization constants Z$, Zp, Z4, 
Zy, as well as the ratios Zp/Zs and Zy jZ A - 

3.2 Outline of analysis procedure 

Our analysis procedure for all ensembles is shortly outlined below. In the 
next subsection, we provide the details for each step and illustrate it with 
an example of the extraction of Zp for ensemble E17.32, see table [2j 

1. Choose ensemble and the correlation function of interest. 

2. For each configuration, average over sites that are equivalent with re- 
spect to the hypercubic symmetry (taking the anisotropy of the lattice 
into account). 

3. For each point X = (x,y,z,t) = (xi, £2, X3, X4), take the ensemble 
average. 

4. Correct for tree-level discretization effects. 

5. Apply some cuts to eliminate points with large remaining discretiza- 
tion effects. This amounts to choosing a direction in coordinate space, 
defined by an angle 9(X) between the position vector of point X and 
the direction (1,1,1,1). The cut consists in selecting only points with 
9{X) < # max > where # max is chosen empirically. 

6. Average over neighbouring directions with equal X 2 = 

7. Compute the Xq dependence of RCs in the X-space scheme Z^{Xq] a, //), 
applying the renormalization condition ([6]) at different renormalization 
scales Xq. 

8. Choose the interval to extract the final value of RC: a <C Xq <C Aqq D . 

9. (comparison with other methods) Convert to another renormalization 
scheme of interest, e.g. the MS scheme, evolve to the reference scale 
of interest, e.g. fi = 2 GeV. 
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3.3 Detailed example — Zp 

We now give a detailed example of the complete analysis procedure drafted 
in the previous subsection. For the convenience of the Reader, we explicitly 
mark the stage of the analysis procedure a given part of the text corresponds 
to. 

Stage 1. We demonstrate the extraction of the renormalization constant 
Zp for the ensemble E17.32. 

Stage 2. The hypercubic symmetry implies that the correlation func- 
tion computed at points which are equivalent from the point of view of 
this symmetry should take the same value. This is indeed true in the free 
theory, but in the interacting case this symmetry is broken in the process 
of stochastic generation of gauge field configurations by the Hybrid Monte 
Carlo method, resulting in different values of the correlation function for 
points equivalent with respect to the hypercubic symmetry. For our analy- 
sis, such points should be averaged over. Given some point (x,y,z,t), the 
equivalent points are the 6 permutations: (x,y,z,t), (y,z,x,t), (z,x,y,t), 
(x,z,y,t), (y,x,z,t), (z,y,x,t). Note that in order to take into account 
lattice anisotropy (T = 2L), we treat t differently - indeed even in the 
free theory the value of the correlation function computed at (x, y, z, t) and 
e.g. (t,x,y,z) is different (unless z = t). In addition to these six permuta- 
tions, the equivalent points include six permutations of all of the following 
points: (L - x, y, z, t), (x, L — y,z, t), (x, y,L - z, t), (L - x,L - y, z, t), 
(L — x,y,L — z, t), (x, L — y,L — z, t), [L — x, L — y, L — z, t). Finally, in 
all of the above points one can replace t —> T — t. This gives in general 
96 equivalent points, although this number may be reduced for points with 
some of the coordinates zero or equal to each other. 

Stage 3. The above averaging is performed for every gauge field configu- 
ration separately. As a result, the correlation functions are fully symmetric 
with respect to the hypercubic symmetry, as in the free theory case. We 
then compute the ensemble average of the correlators averaged over equiv- 
alent sites, together with the corresponding statistical errors. The resulting 
values of the PP correlator (ensemble E17.32) are plotted in Fig. QJleft) as 
a function of X 2 . Even after exploiting the permutation symmetry, there 
still remain points with the same value of X 2 , e.g. (2,2,2,2) and (4,0,0,0), 
leading to additional values of the correlation functions for the same X 2 . A 
possible treatment would be to average over such points and the result of 
such averaging is shown in Fig. [TJJright). The obtained correlator shows very 
irregular behaviour, especially for short distances that are crucial from the 
viewpoint of extracting RCs and making connection to perturbation theory. 

Stage 4. Clearly, the non-smoothness of the correlator at this stage 
comes from discretization effects. Some part of these effects appears already 
in the free-field theory. This is illustrated in the left panel of Fig. [2j Com- 
paring it to Fig. deleft), one can notice the same kinds of structures ("tails") 
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free continuum theory 
M.35, L/a-32, ay=0.00175 i — ■ — i 
P + P + correlator 
non-corr., all points 



(X/ar 



free continuum theory 
P-4.35, L/a-32, afi-0.00175 i — ■ — ' 
P + P + correlator 
non-corr., all points 



IX'ar 



Figure 1: The PP correlation function for ensemble E17.32, (left) plotted 
vs. X 2 averaged over equivalent point obtained by permutations (Stage 2), 
(right) averaged over all points with the same X 2 . 



in the free and interacting correlators. This suggests that it is possible to 
reduce the cut-off effects in the interacting theory correlator by subtracting 
tree-level discretization effects. Following Ref. [54] . we computed for each 
correlator function type Crr a correction factor Ap, defined as the ratio 
of the free correlator on the lattice over the continuum one calculated in 
infinite volume and in the chiral limit: 



A r 



(OrPQOr(0))£r 
(O r (X)O r (0))^ t 



<0r(X)0 r (O))[- c 



7T 4 (X 2 ) 3 



(30) 



where c = 3 for PP and SS correlators, c = 6 for VV and AA correlators. 
Thus, we obtained the corrected correlation functions, C^ r (X): 



Cf r (X) 



Crr(X) 
A r (X) 



(31) 



As an effect of applying this correction to the PP correlator, we ob- 
tain much reduced spread in the values of correlators corresponding to the 
same value of X 2 . We also show the effect of the free theory correction for 
the correlator averaged over points with the same X 2 (Fig. fright)) - it is 
much smoother than the one before the free- field correction (Fig. (TJJright)). 
However, the behaviour is still not smooth enough to reliably extract renor- 
malization constants (note the logarithmic scale for the correlator values). 
This results, naturally, from the fact that the subtracted cut-off effects were 
only the non-interacting theory ones. We expect that the remaining cut-off 
effects, 0{a 2 g 2 ) at leading order in lattice perturbation theory, can cause 
substantial cut-off effects for different types of points (different directions de- 
fined by the position vectors) - hence averaging over points with the same 
X 2 is at this stage unjustified. The role of 0{a 2 g 2 ) effects and their potential 
subtraction will be discussed further below. 
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Figure 2: (left) The PP correlation function in the free-field lattice theory, 
L/a = 32, T/a = 64, a[i = 0. The colour scale shows the angle of the 
position vectors of all points with respect to the direction (1,1,1,1) of the 
hypercube. (right) The (lattice) tree-level corrected PP correlation function 
for ensemble E17.32, averaged over points with the same X 2 . 



So far, our analysis procedure closely followed Ref. [54J. The quality 
of our averaged correlators, after the tree-level correction, is similar, but 
somewhat worse, than the one of Ref. [54J . Clearly, further reduction of the 
cut-off effects that remain in our averaged correlators is desirable. 

Stage 5. For this purpose, we employ an analogue of an idea widely used 
in momentum space computations [U7], [BS], [S3], [ZD]- Here we summarize 
its essence. On the lattice, Green functions depend on lattice momenta 
p^ = | sin Developing p 2 = ^^p 2 ^ in terms of a gives: 

p 2 =p 2 - *L p W + — P W _ -?— p ® + . . . , (32) 
F F \T 36(T 2016(T ' y ' 

where the invariants = X^„pJJ. Hence, for momenta with the same p 2 , 
larger artefacts are expected for the ones with larger pM. The underlying 
reason is the reduction of continuum isometry group 0(4) to the discrete 
isometry group H{4) 0. The orbits of 0(4) are labeled by a single invariant 
p 2 , while lattice momenta which belong to the same orbit of 0(4) do not have 
to belong to the same orbit of H(4). Already for p 2 = 4, a single 0(4) orbit 
splits into two H(4) orbits, corresponding to momenta (1,1,1,1) (pW = 4) 
and (2,0,0,0) (pW = 16). The H(4) orbits labeled by the same p 2 and pW 
can further split into orbits with different p^ etc. (e.g. momenta (0,7,7,0) 
and (2,5,8,0) have the same p 2 = 98 and pW = 4802, but = 235298 for 
the former and 278498 for the latter). 

In momentum space, the above argument implies that in order to have 
minimal cut-off effects, one should consider momenta with a possibly "demo- 



4 For this discussion, we ignore the fact that for anisotropic lattices, the group H(A) is 
broken to H(3), which has minor practical implications for our case of interest. 
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cratic" distribution of the total momentum p 2 over different points, since this 
minimizes the value of the p^ invariant. By analogy, one can expect that in 
position space the points least affected by discretization effects are also the 
"democratic" ones, i.e. the ones close to the hypercubic diagonal (1,1,1,1). 
These are the points with the smallest values of j' 4 l = £^ among the 
ones with a given X 2 . 

That this expectation is fulfilled, one can again observe in Fig. [2]deft). 
For each point, one can define an angle between the position vector of point 
X and the direction (1,1,1,1) - let us denote it 9(X). Different values of this 
angle correspond in Fig. [5] to different colours. We observe a very regular 
behaviour - the points that make up the "tails" are the ones with large values 
of 0(X) and the larger this value, the further away from the continuum curve 
a given point lies. We may suspect that such large 0(X) points enhance 
cut-off effects both in the free theory (as shown in the plot) and in the 
interacting theory. In particular, for such points the tree-level correction 
(|3ip is not enough. Therefore, we chose to avoid them, by defining the 
following criterion: keep only points with 6(X) < # max - Naively, the choice 
of 

$max — should be the best one, as it guarantees the smallest value of 
among all points with a given X 2 . However, such choice is obviously too 
restrictive, as it would keep only very few points, none of them or at most 
one satisfying the restriction a <C X$ <C Aq^ d . Moreover, Fig. E^left) shows 
that the analogy with momentum space is not exact - the most "democratic" 
points lie systematically below the free continuum curve and the "optimal" 
value of 6(X) is around 20-30 degrees, i.e. such values of 9{X) lead to free 
lattice correlator values nearest to the free continuum values. Since the 
free lattice computation is performed at finite L/a, the natural question to 
ask at this point is whether this apparently counterintuitive behaviour is 
related to finite size effects. Simulating on larger lattices, we eliminated this 
possibility - finite size effects are negligible if the coordinates x, y, z, t are 
small in relation to L and T. Since we are interested in the small distance 
behaviour, this condition is satisfied. 

In the interacting theory, the situation is, of course, even more complex. 
It is again true that the points constituting the "tails" are the ones with 
largest values of 0(X). We clearly see that these points are responsible for 
the non-smoothness of the correlator in Fig. fright). Hence, we decided 
to choose 6> max = 30 degrees. This value does not have a rigorous justifi- 
cation, but was found to give the best improvement in the smoothness of 
our correlators while still keeping a sufficient number of points for a reliable 
extraction of the RCs. 

Stage 6. What is more, for most values of X 2 of interest to us, the 
criterion 9(X) < # max selects only points being permutations of (x,y,z,t) 
interchanging spatial and temporal indices (e.g. (y, z, t, x)). Such points are 
not strictly equivalent with respect to the hypercubic symmetry, because 
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free continuum theory 
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(X /a) 2 



Figure 3: (left) The corrected PP correlation function for ensemble E17.32. 
Only "democratic" points were used (0(X) < 30 degrees). Averaging over 
points with the same X 2 involves only points interchanging spatial and 
temporal indices, (right) The corresponding values of the renormalization 
constant Zp vs. the renormalization scale X%. 



of the anisotropy of our lattices in the time direction - hence, we did not 
average over them in the initial stage of our procedure. Thus, it allowed 
us to apply a slightly different tree-level correction for such almost equiv- 
alent points. After this correction, the correlator values for such cases are 
compatible within statistical error and can be safely averaged over. With 
the choice ^max = 30, we do not have to average points that correspond to 
different types of points, like e.g. (2,2,2,2) and (4,0,0,0), since our "demo- 
cratic" cut eliminates such necessity. In this way, we don't mix points with 
possibly very different cut-off effects (such points are then non-compatible 
within the statistical error even after the free- field correction). 

The corrected PP correlator after the "democratic" cut is shown in 
Fig. [3]deft). In accordance with our expectations, the correlator curve is 
much smoother than the one after only free-field correction. Obviously, we 
have not eliminated all cut-off effects and therefore some non-smoothness is 
still present. 

Stage 7. Now, in order to extract any renormalization constant Z-p 
(at some renormalization scale Xq), we impose the X-space renormaliza- 
tion condition directly on the corrected correlation function (using solely 
"democratic" points), i.e.: 



(Cff(X)Cff(O))' (gffpQCff(O)) 
(O r (X)O r (0))*Z * 2 =*Z (0r(X)Or(O))f a r a^.v. 

Since Og(X) = Zf (X )O r (X), we have: 



1. (33) 







7 x (x ] _ C rr (X )^ t _ Crr(X )*£ 
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Using every value of X 2 for which the corrected "democratic" correlators 
are available, we vary the renormalization scale Xq and hence obtain the 
Xq dependence of our renormalization constants in the X-space scheme. 
An example for Zp is plotted in Fig. fright). It shows that Zp is scale- 
dependent (as expected) and reveals that indeed there is still some non- 
smoothness in the data (concealed in the correlator plot, because of the 
logarithmic scale), attributed to remaining cut-off effects, 0(a 2 g 2 ) at the 
leading order. 

Stages 8, 9. In principle, the computed values of the X-space renor- 
malization constants Z* can be used to renormalize the corresponding op- 
erators. In practice, one usually wants operators renormalized in the MS 
scheme, at some reference scale, often chosen to be [i = 2 GeV. Therefore, 
we convert our X-space renormalization constants to the MS scheme, using 
4-loop conversion formulae, as described in Sec. 12.4.21 

For renormalization constants that are scale-independent in the MS 
scheme (Za, Zy), the final values in this scheme are obtained in the fol- 
lowing way: 

^(1/X ) = l/Xo) = Z^ = 2e-^/X ) = Z^(2GeV). 

Eq. <[28j) 

(35) 

The value in the X-space scheme at the scale \i = 1/Xq is converted to 
the intermediate scheme MS at the same scale. Since the value in the MS 
scheme is equal to the value in the MS scheme at the scale \x = 2e _7B /Xq 
and is moreover scale-independent, the perturbative conversion is the only 
step needed to obtain Z v l A (2GeV). Note also that the scale independence 
of Zy t A in the MS scheme and the scale dependence of the conversion factor 
(rather mild, if the renormalization scale Xq is not too small) imply that 
Zy t A computed in the X-space scheme should show an Xo-dependence. This 
is a direct result of the breaking of chiral Ward identities. In practice, it 
turns out that this Xo-dependence is rather mild, as we show in this work 
and was also observed in Ref. |54| . 

In the case of scale-dependent renormalization constants (Zp, Z$) the 
procedure involves additionally the evolution to the reference scale: 

z x {1/Xq) z MS {fi = 1/Xq) = z m^ = 2e-^/X ) 

Eq. H71 

>Z^|(2GeV). (36) 



Eq. J29} 



The value in the X-space scheme at the scale \i = 1/Xq is converted to 
the MS scheme at the same scale, equal to the MS scheme value at the 
scale \x = 2e~ lE / Xq. Using renormalization group equations (Sec. I2.4.3|) . 
we finally obtain z|f(2GeV). 
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p 


a' 1 
[MeV] 


{Xo/af 


Extraction 
£ = 1/X [MeV] 


range 

p, = 2e-~< E 


/X [MeV] 


3.9 


2500 


[7, 13] 


[695, 945] 


[775, 


1060] 


4.05 


3130 


[7, 15] 


[810, 1185] 


[905, 


1325] 


4.2 


3870 


[10, 18] 


[910, 1225] 


[1020, 


1370] 


4.35 


4700 


[10, 21] 


[1025, 1490] 


[1150. 


1670] 



Table 3: Extraction windows for RCs, satisfying the constraint a <C Xq -C 
Aq£ D or, equivalently 1/tt » /J » Aqcd with Aqcd typically about 300 
MeV. In practice, we take Xq > 3a. The values of RCs injiie X-space 
scheme at the scale 1/Xq are perturbatively converted to the MS scheme at 
the scale denoted by fi and are equal to the MS scheme values at the scale 
denoted by ft. 

In order to obtain reliable values of the renormalization constants Z^ s 
at the scale of 2 GeV, the perturbative formulae have to be applied for 
scales that satisfy the condition 1/Xq > Aqcd- Moreover, scales close to 
the lattice cut-off have to be avoided. Together, this yields the condition 
a <C Xq -C Aqcd' which imposes applicability constraints on the method 
- the condition can be satisfied only on lattices with fine enough lattice 
spacing. 

For our ensembles of interest, we have chosen the RCs extraction win- 
dows listed in Tab. [3j Our choices depend on the values of the lattice 
spacing corresponding to each value of the inverse bare coupling /3. Besides 
the physical conditions for which the window has been chosen, the inclusion 
of several points allows for an important check of the scale dependence of 
RCs. The values in the X-space scheme depend on the scale. However, after 
conversion to the MS scheme and evolution to 2 GeV, all RCs should be 
independent of the initial scale up to statistical accuracy, non-subtracted 
O(o?g 2 ) and higher order discretization effects. If this comes out to be true 
in the final analysis, it means that the lattice simulation reproduces correctly 
the quark mass anomalous dimension (Zp, Z$) or the scale independence 
of Zy, Za- Moreover, it provides a strong hint that the extraction window 
was chosen appropriately. For these reasons, we decided to go rather low in 
1/Xq, especially for the case of /3 = 3.9. Whereas 775 MeV might seem dan- 
gerously close to Aqcd , the inclusion of points near this scale does not spoil 
the constant behaviour of MS (2 GeV) RCs converted to the MS-scheme and 
evolved from X-space scheme values at different scales 1 / Xq to the reference 
scale of 2 GeV in the MS-scheme, as we will further discuss below. 

In Tab. HJ we give an explicit example of the conversion/evolution proce- 
dure for Zp. The third and fifth columns show the scale dependence of Zp 
in the X-space scheme and in the MS scheme, whereas the last column shows 



17 



(x / a y 


i/x 

[MeV] 


Zp(l/X ) 


P- 

[MeV] 




r-. ?T\/TQ / ' — N ~r t\ 

Z P ib (2GeV) 


10 


1490 


0.4850(8) 


1670 


0.4685(8) (12) (17) 


0.4921(9) (49) (2) 


12 


1360 


0.4911(6) 


1520 


0.4724(6) (13) (19) 


0.5101 (7) (55) (7) 


13 


1300 


0.4697(10) 


1460 


0.4508(10)(13)(20) 


0.4932(11)(55)(11) 


15 


1210 


0.4734(8) 


1360 


0.4526(8)(14)(22) 


0.5075(9) (62) (20) 


16 


1170 


0.4549(17) 


1320 


0.4340(16)(14)(22) 


0.4924(18)(62)(24) 


18 


1110 


0.4574(10) 


1240 


0.4347(9) (15) (23) 


0.5047(11) (69) (35) 


20 


1050 


0.4589(13) 


1180 


0.4346(12)(16)(24) 


0.5156(14) (76) (47) 


21 


1025 


0.4425(14) 


1150 


0.4183(13) (16) (24) 


0.5017(16)(77)(51) 



Table 4: Example of the conversion of X-space RCs values to the MS scheme 
and their subsequent evolution to the reference scale of 2 GeV. The scale 
[l = 2e~ 7B / Xq. The error for Zp is statistical. Errors for Zp S are: 1st 

- statistical, 2nd - resulting from the uncertainty in lattice spacing, 3rd - 
resulting from the uncertainty of Ajyjg N 2 . See text for more details. The 
chosen ensemble is E17.32. 

that the values converted from different scales 1/Xq lead to similar values 
of Z)F(2GeV), with some fluctuations around the central value. Tab. 2] 
also shows how various types of errors enter at the stage of conversion and 
evolution. The first error is always statistical (for the X-space value it is 
the only error). The second error (systematic) results from the fact that 
uncertainties in the value of the lattice spacing (listed in Tab. [2]) lead to un- 
certainties in the values of the scale 1/Xq in physical units, which is, in turn, 
needed to calculate the value of the strong coupling constant that enters the 
perturbative conversion or evolution formulae. The third error (systematic) 
originates from the fact that we do not precisely know the value of the scale 
Aqcd- To emphasize that our present case is the Nf = 2 one and we work 
in the MS scheme, we denote it by A^jg N 2 - ^ s vame nas recently been 
computed by the ALPHA Collaboration - 310(20) MeV [71] and by ETMC 

- 315(30) MeV [66]. 

Our final values for Zp s (2 GeV) are summarized in Fig. HJlower right). 
The plot shows that within our extraction window, the final value of Z p ls (2 GeV) 
does not depend on the choice of the renormalization scale 1/Xq for the X- 
space scheme. However, fluctuations around the central value are well visi- 
ble. As stated above, the main source of these fluctuations are the remaining 
non-subtracted cut-off effects. We will discuss their role below. 

Fig. flower right) also shows the extraction of the final value of Zp s (2 GeV), 
resulting from averaging over different renormalization points in the X-space 
scheme. To take the systematic errors appropriately into account, we per- 
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Figure 4: The values of RC Zp in the MS scheme at the reference scale of 2 
GeV. Ensemble: E17.32. We show the influence of applying tree-level cor- 
rections and "democratic" cuts on the final values of Zp S (2GeV) obtained 
from the X-space scheme at different scales Xq. Upper plots: no tree- level 
correction applied, lower plots: with tree-level correction. Left plots: no 
"democratic" cuts (all points), right plots: only "democratic" points. The 
exclusion of "non-democratic" points from the averages results in vertical 
shifts of some of the points; hence, the data points shown in the right plots 
are not merely subsets of the points shown in the left plots. The final plot 
of our analysis for Zp for this ensemble is the lower right one. Errors as 
explained in the caption of Tab. O Note that the vertical and horizontal 
scales are the same in all plots. 

form this averaging in the following way. We take the central values of 
Z^ s (2GeV)(l/X ) (where this notation shows that they were extracted 
from different renormalization points in the X-space scheme) with their sta- 
tistical errors and compute the weighted average. The error of this weighted 
average is then the statistical error of the final value of Zp S (2GeV). Next, 
we compute the weighted averages for the values of Z P iS (2 GeV)(l/X ) 
shifted upwards and downwards by the lattice spacing systematic error. The 
spread of the resulting weighted averages defines the lattice spacing system- 
atic error of the final value Zp IS (2GeV). Finally, we similarly extract the 
Ajyfg N =2 related systematic error. The ultimate error of Zp s (2 GeV) comes 
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from the statistical and two systematic errors combined in quadrature. Such 
procedure finally gives ZjF(/* = 2GeV) = 0.5033(10) (59) (17) = 0.503(6) 
for ensemble E17.32. 

The different plots of Fig. 0] illustrate the role of tree-level corrections 
and "democratic" cuts. The upper left plot corresponds to non-corrected 
correlators and using all points. Depending on the choice of Xq, the X-space 
RC value converted to the MS scheme and evolved to 2 GeV can yield values 
that are up to ca. 25% away. For coarser values of the lattice spacing, this 
difference can even exceed 50%. Therefore, as we have argued above, the 
tree-level correction and "democratic" cuts are crucial for a reliable analysis. 
As the lower left and upper right plots of Fig. [5] show, the correction or 
cuts exclusively provide a big improvement over the original correlator case. 
However, it is obvious that only the two corrections combined can give a 
really good quality of the extracted value of Zp. 

This concludes our detailed discussion of the extraction of Zp for ensem- 
ble E17.32. In the next section, we summarize our results for all RCs and 
all ensembles. 

4 Results and comparison 

4.1 RCs for ETMC ensembles 

The procedure for other renormalization constants follows the line of the 
detailed example of the previous section. To compare the quality ofthe data 
for different types of RCs, we plot in Fig. [5] the dependence of Z^ s (2GeV) 
on the X-space scheme renormalization point. We show the results for Zp, 
Zg, Zy, Za and the ratios Zp/Zs, Zy /Za for ensemble E17.32. The errors 
for Zp and Zg are much larger than the ones for Zy and Za-, since the 
latter are scale-independent in the MS scheme and hence do not require the 
step of evolution to the scale 2 GeV. For the ratios Zp/Zs and Zy/ZA, 
we quote only the statistical error, because both the RC in the numerator 
and the denominator have the same conversion factor and the same scale 
dependence - hence, the conversion factors cancel each other and similarly 
for the scale dependence. Thus, Zp/Zs and Zy/ZA are fully scale- and 
scheme-indep endenlQ. 

We now move on to discuss our results for other ETMC ensembles of 
gauge field configurations. In particular, for two values of the lattice spac- 
ing, we can perform an extrapolation of RCs to the chiral limit and obtain 
their mass-independent values. Such extrapolations for f3 = 3.9 and 4.05 
ensembles are shown in Fig. [6j In all cases, our results are compatible with 
a linear dependence on the quark mass and this dependence is rather mild 

5 The RCs Zv and Za are often said to be scale- and scheme-independent. However, 
this concerns only renormalization schemes that preserve chiral Ward identities. 
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Figure 5: Extraction of final values of RCs Zp, Zs, Zy, Za and the ratios 
Zp/Zs, Zv/Za for ensemble E17.32. Errors as explained in the caption of 
Tab. HI Note that the horizontal scale is the same in all plots. The vertical 
scale is always [Z-p — 0.075, Zp + 0.075], where Zp is the extracted value of 
RC. 

(in some cases the slope is compatible with zero). We also observe that the 
chiral limit value is always compatible with the value at the lowest quark 
mass. Therefore, at /3 = 4.2 and f3 = 4.35, where only one ensemble is 
available to extract X-space RCs, it is plausible to use our extracted values 
as mass- independent values, since they correspond to the same pion mass 
as the lowest quark masses at j3 = 3.9 and 4.05. 

Tab. [5] summarizes all our final results for the renormalization constants 
extracted in the X-space renormalization scheme, together with the values 
determined in the RI-MOM scheme, given for comparison [22], [34J. Al- 
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Figure 6: Chiral extrapolations of renormalization constants for f3 = 3.9 
(left) and /3 = 4.05 (right). 



z P 



z s 
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Zp/Zs Zv/Za 



B40.24 
B64.24 
B85.24 
B150.24 
chiral 
[22] 
3-1 



0.458(2)(12)(21) 
0.462(2)(13)(23) 
0.463(2)(12)(22) 
0.460(2) (12) (21) 
0.460(2)(14)(25) 
0.437(7) 
0.457(10)(16) 



0.765(3)(21)(38) 
0.776(3)(22)(41) 
0.778(4)(22)(41) 
0.807(7)(24)(45) 
0.750(5)(25)(47) 
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0.626(1) (2) (6) 
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0.872(5) 
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C30.32 
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0.695(10)(13) 


0.702(6) (1)(1) 
0.686(1)(1) 


0.785(7) (1) (2) 
0.789(1)(2) 


0.707(14) 
0.713(2)(2) 


0.893(16) 
0.869(2)(4) 


E17.32 


0.503(1)(6)(2) 


0.681(2)(8)(3) 


0.692(1) (1)(1) 


0.768(1) (1)(1) 


0.740(3) 


0.902(3) 



Table 5: Renormalization constants for all ensembles, extracted in the X- 
space renormalization scheme, converted to the MS scheme and evolved to 
the reference scale of 2 GeV for scale-dependent RCs. For (3 = 3.9 and 4.05, 
we show our values in the chiral limit. For comparison, we also give values 
determined in the RI-MOM scheme [22], [34j . Errors are: 1st - statistical, 
2nd - resulting from the uncertainty in lattice spacing, 3rd - resulting from 
the uncertainty of A^jg N =2 - The last two errors do not apply for Zp/Zs 
and Zv/Za, see text. 

though the cut-off effects in the X-space method and in RI-MOM are, in 
principle, very different, in most cases we find compatible values. This 
strengthens the expectation that both methods should lead to only rather 
mild discretization effects. The errors of the X-space method are in gen- 
eral slightly larger than the ones of RI-MOM. This is particularly visible 
at (3 = 3.9, where the by far dominating error is related to the necessity 
of going rather close to A^jg N _ 2 (down to around 800 MeV), where the 
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B40.24 
B64.24 
B85.24 
B150.24 
chiral 



0.458(24)(25) 
0.462(26) (26) 
0.463(25)(27) 
0.460(25)(26) 
0.460(29) (29) 



Z S 

0.765(44)(73) 
0.776(47)(78) 
0.778(47) (76) 
0.807(51)(90) 
0.750(53)(90) 



Zv 

0.626(6 
0.629(6 
0.629(6 
0.648(7 
0.616(8 



(31) 
(33) 
(34) 
(43) 
(40) 



Z A 

0.729(8)(15) 
0.725(8)(14) 
0.725(8)(14) 
0.725(8)(14) 
0.729(9)(17) 



Zp/Zs 
0.599(4)(35 
0.595(5)(36 
0.597(5)(33 
0.574(8)(41 
0.609(6)(42 



Zy /Za 
0.861(4)(25) 
0.869(4)(28) 
0.872(5)(33) 
0.900(9)(41) 
0.847(6) (34) 



C30.32 
C60.32 
C80.32 
chiral 



0.469(9)(13) 
0.475(10)(15) 
0.480(10)(18) 
0.463(16)(23) 



0.705(16)(40) 
0.709(17)(39) 
0.724(18)(42) 
0.692(27)(68) 



0.647(3 
0.648(3 
0.652(3 
0.644(4 



(19) 
(20) 
(23) 
(34) 



0.731(3)(12) 
0.732(3)(11) 
0.743(4)(11) 
0.724(5)(19) 



0.669(5)(20 
0.672(6)(15 
0.668(6)(14 
0.671(9)(31 



0.886(4)(17) 
0.886(5)(15) 
0.879(5)(19) 
0.890(8)(29) 



D20.48 0.516(8)(10) 0.730(11)(17) 0.702(6 



(17) 0.785(7)(10) 0.707(14)(8 



0.893(16)(13) 



E17.32 0.503(6)(12) 0.681(9)(16) 0.692(1 



(15) 0.768(2)(11) 0.740(3)(3) 0.902(3)(14) 



Table 6: Comparison of the total error for each RC extracted in X-space 
(the number in 1st parentheses) and the magnitude of fluctuations of RCs in 
the extraction window (2nd parentheses), due to the un-subtracted cut-off 
effects (0(a 2 g 2 ) at leading order). The total error is the statistical error and 
two systematic errors, combined in quadrature. 

strong coupling constant and hence the conversion factors X-space— > MS 
strongly depend on the actual value of N f =2- This uncertainty becomes 
less important at finer lattice spacings and then in some cases we can achieve 
comparable or even smaller errors than RI-MOM (e.g. Zp, Z$ at (3 = 4.2). 
In general, our precision in Zp, Z$ increases from 5-7% at the coarsest 
lattice spacing to better than 1% at its smallest value. For Zy, Za it is 
approximately 1% already at f3 = 3.9 and it reaches even 1-2 per mille at 
(3 = 4.35. We observe a similar behaviour for the ratios Zp/Zs an d Zy /Zj^. 

4.2 0(a 2 g 2 ) effects 

As we have discussed above, our analysis procedure involves the subtraction 
of tree-level discretization effects and "democratic" cuts, which eliminate 
points that are subject to enhanced discretization effects due to the break- 
ing of rotational symmetry. However, we have not eliminated any of the 
leading order of lattice perturbation theory 0(a 2 g 2 ) effects. As was already 
suggested in Ref . [5l] , the fluctuations of RCs observed in the extraction win- 
dow a <C Xq <C Aqq D can be due mainly to these cut-off effects (and higher 
order ones). In other words, if true, subtracting these effects as a part of the 
analysis procedure should lead to compatible values of RCs for all choices of 
the renormalization point Xq. In our study, we found further evidence that 
the fluctuations are indeed caused by un-subtracted discretization effects - 
the foremost being that their magnitude clearly diminishes with decreasing 
lattice spacing. Moreover, the effects seem to be rather regular, when one 
compares the same Zy for different values of j3 - suggesting that these are 

6 Note that /3 = 4.2 is an exception to these observations. This is only due to statistical 
accuracy - for the ensemble D20.48 we have only limited statistics. 
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not random fluctuations, but some systematic effects. 

Our estimates of the size of these fluctuations is shown in Tab. [6l which 
shows the extracted values of RCs, together with their total errors (statistical 
+ two systematic errors combined in quadrature) - the number in the first 
parentheses, and one half of the spread between the highest estimate of 
Zp IS (2 GeV) in the extraction window and its lowest estimate - the number 
in the second parentheses. In general, the magnitude of fluctuations is 1-15 
times the total error, with a systematic tendency to decrease for finer lattice 
spacings. 

A systematic way to improve the X-space method would be to compute 
these effects in lattice perturbation theory or numerical stochastic pertur- 
bation theory. This has the potential to reduce the size of fluctuations and 
hence enhance the range of applicability and the reliability of the method. 
Our present approach averages the fluctuations out. However, since, 0{a 2 g 2 ) 
effects can be small for some types of points and much larger for some others, 
it will be important to understand the size of these perturbative corrections 
in detail. 

5 Conclusions 

In this paper, we presented a feasibility study for the application of the fully 
gauge invariant X-space renormalization scheme method to extract renor- 
malization constants of bilinear, flavour non-singlet operators. To this end, 
we have used ETMC Nf = 2 ensembles of maximally twisted mass fermions 
on a tree-level Symanzik improved gauge action. Using continuum pertur- 
bation theory formulae to convert between the X-space scheme and the MS 
scheme, we also computed the values of RCs in the latter, at a reference scale 
of 2 GeV. This allowed us to compare our results to the ones obtained in the 
RI-MOM scheme and likewise converted to MS. In general, we found very 
good agreement between the two non-perturbative schemes, signaling rather 
small cut-off effects in both approaches. For our coarsest lattice spacings, 
our errors are visibly larger than the ones of the RI-MOM scheme. How- 
ever, going to finer lattice spacings, the errors become comparable in size. 
Hence, as expected, we conclude that the X-space method works best for 
fine lattice spacings, where the renormalization window, the region of small 
distances satisfying the constraint a <C Xq <C Aqcd> becomes sufficiently 
wide and allows for a reliable extraction of RCs. In practical terms, the 
lattice spacing of approximately 0.08 fm (our /3 = 3.9 ensembles) is on the 
verge of applicability of the method. 

The X-space method was successfully used before for the computation 
of RCs in the quenched case with the standard Wilson gluonic action [54] . 
while in this paper, we work in a dynamical setup. Apart from this, the main 
differences with respect to Ref. [5l] are: the use of "democratic" points to 
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reduce cut-off effects and 4-loop conversion formulae between the X-space 
and MS schemes. Both seem to be important ingredients for the reliability 
of the X-space method. Further improvement can presumably be obtained 
by calculating 0(a 2 g 2 ) effects in LPT or NSPT and subtracting them from 
the correlators. Thus, the values of RCs would be even less contaminated 
by cut-off effects. 

To finalize, let us summarize the main features of the X-space method. 
Being fully gauge invariant and free of contact terms, it prevents mixing 
with certain types of operators. This makes it suitable for the computation 
of RCs needed to renormalize weak matrix elements. This, however, requires 
the additional step of deriving continuum perturbation theory expressions 
relating the X-space RCs to the MS ones. The X-space method has also 
some practically appealing features - it is easy to implement and relatively 
cheap computationally. Thus, we believe that the X-space renormalization 
scheme has promising potential as a method to perform non-perturbative 
renormalization in lattice QCD in the future. 
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